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The  basic  equations  required  to  acquire  an 
understanding  of  the  autopilot,  roll  system  and 
guidance  computer  for  a  homing  missile  are  derived. 

* 

H cm  the  use  of  these  equations  lead  to  methods  for 
analysis  and  design  of  the  systems  is  shown. 

The  factors  that  affect  the  desif.ni  are  enumerated 

and  the  most  important  ones,  explained  in  detail. 
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lo  INTRODUCTION: 

This  report  has  been  written  to  serve  possibly  two  purposes,  First  to  arn»a-ir,t 
people  not  familiar  with  the  autopilot,  roll  system  or  hcminp  puidance  system  vlth 
IT^°1  Pr0b\T  inVOlVed  in  their  desiE"  *"d  suggeSnSiSs  of %£££  for 

I  ork^Lth1seh^fr’lemSL an?  s?cond„to  presenfc  under  one  cover  some  of  the  analytical 
bE^i0  “t*®  desien  bf  the  systems  o  There  are  therefore  both  elementary 

No  menSof is^S^1011  regafllns  the  s^ster,s  and  methods  for  analysis  and  synthesis. 

fo-  IJemne  The  iLZanS;0?UG^-mp,iter  techn^ues  sl"ce  «*■  i«  familiar^round 

Sots*  shoula  •" usem  ~ 
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2o  AUTOPILOT  DESIGN: 

The  autopilot  will  be  defined  as  the  comole to  system  for  which  the  input  and 
the  output  are  the  following: 

Input  -  Coraaand  acceleration  in  any  directi on  normal  to  the  missile  center** 
line. 

Output  «*  Resultant  acceleration  of  tho  missile 0 

The  airframe  characteristics#  the  control  systems,  and  the  sensing  instruments, 
if  any,  are  involved. 

2<>1  Derivation  of  Transfer  Functions  For  The  Airframeo 

The  basic  design  of  the  autopilot  cm  be  carried  on  by  considering  tho  system 
to  have  three  degrees  of  freedom  (l)  missile  rotation  in  pitch,  (2)  missile  Cogo 
translation  normal  to  the  body  in  the  pitch  plane  and  (3)  control  surface  ds flection 
to  produce  motion  in  pitch  plane. 

Two  additional  degrees  of  freedom,  (l)  rotation  of  missile  about  longitudinal 
axis  and  (2)  tail  surface  deflection  to  produce  roll,  are  included  in  the  final 
analysis  to  chock  the  compatibility  of  the  autopilot  design  with  the  roll  control 
system  in  the  presence  of  aerodynamic  coupling.  The  equations  of  motion  and  the 
transfer  function  for  this  analysis  are  developed  in  the  section  on  roll  system  design# 

The  angles  for  the  single  plane  analysis  are  shorn  in  Figure  lo  The  sign 
convention  is  also  shown. 

The  Angles  are  the  Followings 

°<  “  Angle  of  attack  *  angle  between  missile  centerline  and  velocity  vac tor 0 
“  Angle  between  the  velocity  vector  and  the  reference. 

“  Angle  uetween  the  missile  centerline  and  the  reference, 

1  *  Angle  between  the  control  surface  and  the  missile  centerline o 

The  relationship  between  the  rate  of  rotation  of  the  velocity  vector  and  the 
acceleration  of  the  e.go  normal  to  the  missile  body  v&ll  be  derived  first. 

“A* 

Let  j  be  a  unit  vector  along  V  and  let  K  be  a  unit  vector  -L  to  ^  and 
positive  in  the  ^direction  for  increasing  If  »  (See  Figure  l)9  then  the  velocity 
of  the  Cog.  is  V'V  X  and  the  acceleration  of  the  e.go  is 


i 

i 


dl  «  V?  -h  v  —  bor,  dZ  m  Y  it  (y  w  R«#/$trc) 
dft  dt  c/t 


3?1JL_A_JL 

"r-yy  C& 
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therefore, 

Soi-l  A]/  -  '/J  -f  l/J  ?  (^AeA 

PC 

V  '$  is  therefore  the  acceleration  of  c„go  normal  to  tne  velocity  vector 0 

The  acceleration  of  the  c&g»  normal  to  the  missile  body  is  defined  by  H? 

(g’s)  anu  along  the  missile  axis,  by  ($'*)  •  ^ 

The  relationship  between  the  accelerations  with  respect  t  the  missile  axis  ana 

•  f 

^  iSz  ~  ^3  op 

e  » 

if  ^  is  in  deg/sec  and  fl#  oC  is  assumed  to  be  negligible  compare^  to  On  CrSi-  °v 


14  ~ 

WV-sr  }j 


i2*r». 


This  relationship  is  used  frequently  in  l he  following  sections,, 


The  aerodynamic  forces  normal  to  the  missile  body  /?*  is  a  function  of  c<  arri  l_ 
r  0 
\  J  /  .m  _ .  /"  /  *  /  I  t  /  __  .  I.  •/  <*t  /  if .  /y  jii  t  b  .4 


W  fa  T  (**}  0  where  W  -•  usdCfh  {  d  **  pQunl  $ 

o 

\J/  J/t7  C;  J_f  A/  J  1  S’  Y  \ 


Assuming  the  partials  to  be  constants  and  integrating 

*  VI/  /?<?  «  ^  <Lf  ^ 


2, 1-3  VI/  /?*  -  ol/  ^  4-  J  f  i 

Similarly  for  the  aerodynamic  moments,, 

a.i-v  J  £  -  4**  -f 

O  ^  C/  * 

X  *  pitch  moment  of  inertia  (slug  fto^).. 

i 

The  moment  is  also  a  function  of  ^  but  this  term  is  neglected  since  the  control 
system  of  the  type  proposed  wifco  rate  gyro  feedback  makes  it  :.  negligible  factor < 

The  partial  derivatives  4~L  ^  ^  and  -LUtt  are  one  fern  of  the 

stability  derivatives,  Normally  toey  are  written  in  the  fom 
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-ft  g  (J£s) 

Ch'=  t*  it  (%) 

c-'  &  7^(  y^j) 
tm^m&  ^C'/u) 


*»  dynamic  pressure  «  l/2  ^  ( slug/sec d  ft«) 

~  }*+?l  Pi  A/  * 

where  *  density  of  air  (slug/ ft"?) 

PK  *  ratio  of  static  pressure  at  altitude  to  that  at  sea  level* 
ft \  “  Mach  Poo 

5  «  reference  area  (normally  body  area)  (ft-2)  “  «► 99li  ft.2  for  Tartar 

d  "  reference  moment  arm*" normally  body  diameter)  (ft*)  «  1*125  ft.  for  Tartar 

For  purposes  of  analysis  the  force  and  moment  equations  2*1-3  and  2oI~ii  are 
put  in  the  form 


2„1«5 


cf  'ss  A*  +  Si 
^  »  Coc  -b£  l 


where  all  angles  are  in  degrees. 

The  relati  onsaip  between  taese  coefficients  and  the  original  stability 
derivatives  can  he  easily  ded  ced„ 

A  A 'isr+rVwg’/)JU  M  *G>- 

,  vw 

b  *  Ctt . 

\/lV 

c »  Cl zjjtiyy/  Ll m3( 


•r*** 

8  “  (*Jj3u)lum\3LLcL  N  Y'm 

X 


4 


1  • 
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The'  values  of  C  Wt<:  and  must  be  the  'ialues  for  the  particular  Ccg» 

location  considered..  If  these  coefficients  are  available  for  a  reference  c0g0 
location  and  have  the  values  0Kif  a,tej  (%t  ,  the  values  at  any  other  c<,g„  location 


are  given  by0 


c. 

0 


-  C' 


-  ''  m 


®C* 


-£  X"—  /o  si 

~7T 


m *  ~  t'  -t 


V5^  ^ni 


where  XQ  is  the  reference  cego  and  X  is  the  desired  c0go  Both  these  quantities  are 
in  inches  and  measured  front  Station  0  which  is  at  or  near  the  nose  of  the  missile „ 
d  is  the  reference  moment  arm  in  inches  (l3e5.  inches  for  the  Tartar) „ 

The  transfer  functions  for  the  airframe  can  be  derived  from  equation  20l-3° 
Assuming  all  initial  conditions  are  zero,  the  operational  form  of  these  equations  arec 


Co 


1«6 


"5  <r  ~  fy  ~t  B  2 
~  Cx.  +  £  i 

from  the  definition  of  the  angles 


substituting 


6r 


<5 

5  v  -  /I  -t  /}  -=  B  l 

-CV  -f  O  «  £  ? 

(  ■$  I1- ft)  y  /i  'if,  rs.  8  1 

C  flr  -t(s--c  ]  tJ,  =  fz  t 


By  use  of  Cramers  iiule 


lB  -fi  j 

I  £  (s  l-c  ) 


C  htc) 
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2„1«7 


2.1-8 


— -  3<C  ~t  Alz~ 
ft  $  —  c. 

-  d£d?£'  \ 

'  c-T“  )  ~~ 


V  tr  —  /  /?C  \ 

^  i  “cT  ; 

/5^  Bl 

'i  •*  /  C  erl 


rter-Si 


a  ^  5  +  - 


s*h  4  s  ~c 

-  JL^jLA?  -PC. 

~r 

Ah  . 

%  -*  ~/fl£-BC\  l  -+■  -----  S' 

'  l  (.  ~ 7T- j _ />e-gc 

v-  -_sv  _1  ;r;  :r  ■■/• — 5 


^  ~  -VV  -~/uTs=7-  / 

Substituting  /'/  for  y-  equation  2.1*7 
0  5 

?ol~9  n^A  «  — dr 


ypvi-  '  c  y 

For  the  Tai*tar  missile  the  values  of  A  anc  C  -//ill  vary  with  angle  of  attache 
To  facilitate  the  analysis,  however,  constant  values  for  these  coefficients  are  used. 
The  analysis  will  be  valia  for  the  angle  of  attack  region  for  which  the  coefficients 
were  picked. 

The  values  for  these  coefficients  for  the  Tartar  missile  are  t  ■.  *en  ir  the 
following  table  for  two  typical  flight  conditions* 


Case  Number 

Flight  Conditior 

&<1 

A _ 

B 

_ C  _ _ 

_ B _ 

1 

r  k  3.s\ 

6° 

lo7S 

MS 

*>32 

'•391 

•< 

2 

Isa  ) 

18° 

2- .2 

MS 

oO 

-391 

3  ( 

"n  s  ,v 

6° 

*6$ 

c-  226 

-X 

Ml 

30,000  fto  f 

Altitude  J 

10° 

lo0£ 

,226 

b 

Ml 

These  valves  are  for  a  corbined  plane  maneu.er,  i„e.,  the  resultant 
angle  of  attack  and  g’s  are  in  a  plane  1*5°  from  the  plane  of  the  i  ingsc 

Substituting  the  values  from  the  table  in  equations  2.1-6  and  2, 1«»9 
these  eouations  became,. 
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{ .f  < 

4y  *  6  - 

.7  *■>-•  J 1  ;\  O^'-d '  -J-  f 

lH  '"sr  '■Z&kj'  - .  /  ~  ■-  00S2&SSt~ 

tS^/sr  .  ,**2  j  S‘ , 


W  _  ou,  ._.„  ;  -/•  ,  ^ _ 

/i  ~  ■  -  ,  OJ  7X0  V-  .’ 

#*/  y-  <?fV/  /  .  o 

~v ohj-^  as.vs'S-o  j 

•'  _ _  /  "T"  f.tyi'J? 

A  y  — 

'Vs  "ft-  Iiy  - s  *• _ 

*i  /#v5"  ~  •  1*7  jT  „*/  v’  / 

>7  — ~-^T-,&£4. _ _ 

-  -,Zs&  ■ 

■  n%*3£>U  -—■/£■:?■! - 

*’c  'Wvs'"  -a  t 

These  equations  are,  in  general,  of  the  form 
'•/)  =  K  .  /v  7v'-> 


VI  -  * 


/-/'f>i7'5  /-  T3 


‘>1/ 

'/  •«  J£_  V  I 

*  0  J<?/Lc~  <'V  — . ~~~ 

/tZi&m* 

• ‘/ 

)  ms-  ( • 1  *V  =v 


The  significant  factors  that  are  apparent  from  these  equation?  are: 

lo  The  airframe  is  essentially  s  second  order  syetera* 

2o  The  airframe  by  itself  can  be  di\ir gently  unstable  for  certain  flig.it 
conditions  and  angles  of  attack.  For  cases  2  and  4- ths  roots  of  the 
characteristic  equation  have  positive  real  partSo 

3-  The  zero  frequency  gain,  K,  varies  an  magnitude  and  polarity  v.itn  flignt 
condition  and  angle  of  attach 

I 

Liv  is  approximat.-ly  the  derivative  of  0^  o  That  is  it  can  be  used  as  a 

measure  of  the  rate  of  <?»  4  ^  feedback  can  therefore  be  used  to  provide 
dsnping  for  the  system,-  ^ 
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2o2  Form  of  Control  Equations 

The  control  equation  can  hevo  any  number  of  forms*  The  most  common  are  the 
following t 


2*2*3. 


K  ii 


(  Tupe  A  \ 


2"S,'J  U  -  -Kile.  rn,^  s) 

2*2-3  lt  *  |C,  (n  -«<.•)  4  If*  «/'  ("/yr  c  } 


2c2-U 


ic  -  j  '.IhJt :  (Typ*  />) 


where  4.  is  the  command  wing  def 2  ection  and  /?c  is  the  command  acceleration  jand  the 
prime  represents  the  instrument  outputs* 

Block  diagrams  of  the  complete  autopilot,  sys  tem  using  these  control  equations 
are  shorn  in  Figure  2  «*  (a),  Co),  (c)  and  (d}„  These  systems  wilL  be  explained  in 
detail  after  the  closed  loop  transfer  f arctic  >s  xor  eacn  are  dei*i/ed~ 

The  closed  loop  transfer  fu  'Ctions  car  be  oerived  by  using  tie  equations  from 
section  2-1,-.  I  he  foil  owing  a&ecroptions  >  all  be  irate  fo->'  the  senro  anc.  instrument 
responses* 

CrO)  *  1 
6,\  ($)  ■=  / 


Cu  ft)  *  *  h>-  £ 

“  5  F/t  Type-  D 


These  assumptions  do  not  iu validate  the  general  conclusions  that  can  be  drawn 
for  these  systems*  The  actual  choice  of  gains  however  cannot  be  made  without 
consideration  of  these  transfer  functions* 

For  the  type  A  system,  substitution  of  equation  2„1~9  into  equation  202"1 
results  in 


2 -2-5 


V  =  kv  /  ite-ec  \  /  4  _£ 4- 

^  /TTr  <  'T-  J 

-  ^  5  4/ 


X-AJ5S& 
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For  the  'remaining  types  of  control  systems*  aquation  2«1»6  is  used  -vdth  tbs 
control  equations  •. 

For  type  B  the  derivation  is  as  follows.: 

(  '  /  —  ,.y  ~  JH  /  -  t-’1 

C-,y  -y(s'-c)t  ~.£  I  * 

Kz  V  -  t  K  *  ' 


-  a 

*\e. 

-  1 

r> 

-'c 

••  F 

— '  t’t 

_ K,  (  -•  8C  i-S  s1;) 

~  S  6- c  t  -£fl  * 30  -  6  ="  hfi7y-'il) 


2o?™5  ,,  ^  y„.  £ 

I’Aclif  • 


L i _  _ 

t '  **«■ 1  «•'  “/  -f"  -  t  txtjTTn  W(mF-b£) 


For  type  C  toe  derivation  is  as  follows; 
(r>  v-  -  4  V'  -  ,3  3  r  <? 

c  r*  f  i  ■*  o 

o  .y  <■  ^  .t  5/  c  y  „  ?  -  j' . 
%1<S r  1  '  /w  1 


<  j‘  i  ■ 

y  .  ?  * 

!  0 
o 

-0  -B 

j»*i  -  £. 

1  r'« 

KzS  - 1 

<S 

•../}  -  6 

C 

•■  v 

V.F. 

■wvfr* 

f„s  -*■ 

*  V 

2.2 -7 


8  s ----- 

k  V  8  ‘75  _ 2* 


(*&&  -t)s  ;  (/4£r“p)  _  Ki.  C 

fas-M/iB-Bc.)  K/Af-SO  J-K'  jfeK,(t  e-BCi 
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(s  +#)  ir  -  a  y  -  8 1  ~  o 
c  <r  +  ( [s'^-c  )  ->//  -£  \  ~  o 

t  (KlS  S1  ■  K/  >v 


tP¥S" 


0 

- 

-6 

0 

s%-c 

-£ 

3  S  “ 

-s 

1  Sr.*? 

- 

-6 

C 

.«  1  r<: 

-7 

f&V 
'if HS 

S  Kj.  St  ft  s  *• 

~S 

„ .  tfs  *■  ~*r  * 

_ l±Jg[Zg£ZEcl-- 


2,2-3  "A/>\r 


^ . _  1 

kMM-ZQ  k'V  fae-Bt )  &v( rtf-M)  K^t. 

/^v-r-  /res*  I  mr} 


/»¥& 


Fo r  purposes  of  ex&mi  ruling  t.ae  n/  Y\c*  transfer  .functions;  qualitatively 
the  following  appi'oxiEatiuss  will  be  wade 

f\E  -6  C  /9£ 

'prs^irc:  —  c  . 

The-  transfer  functions  for  ine  four  types  of  systems  become* 

Type  A  (  l6  -  “  K,  ) 

2o2"’9  wtfGt2*.u..  /  —  Cl~&,  rf  --~y  • . -—■— — — : »»jt~ - -  — 

/n<-  iiftf-s  \  c.  ?  — -S*  -  .Ely*-  / 


Type  B  ( Tc  *  ~K,  v\c  +  ?C  ^  ) 

2o2“10  V  ^  J£lX  f  -4^- )  --npr — 7--— --/"»-y- - - - - rr 

/r'c  /*»vur  w-  >  ~-y  -  f  M  +  ili  r  -  >  ’  y-  is* 

■fype  C  (ic  *  (V  “  ^  )  +-/fA  f  1  ) 

2*2-11  Wu  :>*  „ - t— - 7 - £ - , - - - - - 

r-  -  -S  t,  (thE-fi)  ..  &i 

~  ■*  73'7.r  /T*  /feaT 


-VL,  .v  £  /* 

/^/.r  " 


rvST 


H 


-■  'W:  ft  H  V  A-  1  R 
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2n  24.312  %c  * 


-5 


-  ~> 


-A  S-y- 


‘>u  (l0%  -/  42 \  {  Jz ,  \ 

(sf£r  &^Tr  )s> ('  *  TsTv- } 

•  '***  '*w)  \  tr&J 


&£££.  M£d£ 

/5>'/jr  /Wr  '  '-w  -W**/  l  Twrj 

"  vP®  ^  is  the  simplest  forii  poaslJlea  This  oysttm  can  operate  only  «*er  ih« 
rcnga  of  angle  of -attack  where ^  is  fairly  constant  and  the  sirfcvana  ,is  ©table 

V  feoalfo  Tho  gfin,  %  must  la  matte  fco  vary  with  flight  conditions  ii  the  aero 
irequency  gain  (static  gain)  is  required  to  ha  constant*  The  damping  £«•  the 
astern  is  provided  only  from  the  airframe*  l»;  general  the  damping*  odef fieiunt  P’ 

,  ii»aoWoa  of  critical  damping)  is  on  the  odder  t£  *1  to  *2  for  thosu  angles  of  tuack 

\vnera  ihajzJL.iiS  approximately  constant* 

,*  '  -o  , ,  ■  .,  "  '  **  — r  * 

Type  B  uses  i\  rate  gyro  feedback  to  improve  ive  daisoing  cViuraov eristic  of  the 
wjB.tma  The  requirements  on  "£-0  stability  cf  the' airframe,  and  K-..  are  the  same 

as  for  the  Typo  A  system  Kg.  can  be’  chosen  to  produce  any  <f  desired, 

Type  C  system  uses  an  additional  accek-cmster  feedbacks  This  reaves  the 
restriction  on  variations  in  and  the  stability  of  the.  atsSraais  *£  .the  gains  Ki 

and  Vp  are  chosen,  properly*  If  it  tl  desirable  to  maintain  •  aoorosdnately  the  >ame 
apseo^of  response  and  damping  characteristic  at  *11  flight  conditions,  these  gains 
can  ba  made  to  vary  vithsorr.o  maasarshie  quantity  which  varies  in  tbs  same  manner  ■ 
os  tne  aerodynamic  ftoefflcient»«  Total  pressure  is  such  a  quantity,, 

The  variation  in  static  gain  with  variation  In  <j  will  depend  on  tix«  highest 
vcj.ua  vnat  can  be  used  for  Kb0  Ths  restriction,  on  the  highest  value  for  is 

explained  in  die  rte:ct  section*  •<• 

•  6 

^  s^aton;  in  similar  to  th sit  ised  fox'  the  l’artsz,.  In  this  ovsteiu  a  rats 
feedback  control  surface  servo  is  used*  This  results  in  on  additional  intesratiod 
and  requires  a  feedback  term  porpo-'Ucnol  fco  the  derivative  of  the  rata  gyro  sxgual. 
to*,  autopilot  is  now  a  third  order  system  instead  of  the  .second  order  nvoweflw  for 
m,e  other  three All-three  coefficients  for  the  3rd  order  charcctov.ia-.Lc  cautoiro 
can  bo  oonti‘oll.ed  by  pruplr-chJoic--^l5l  Kp  &Td  K*,  Tho  sfci>  gam  ft*  this  ayutaa 
io  keen  uo  (xepand  only  «  3.f  thiS’  uwr.r.  Ife  apprexima  y  n  constant,  than  abs 

static  gain  win  be  approxife-ly  constant  for  all  angles  of  attack  and  flight 
conditAonRo  If  the  term  .is  not  approximately  c  constant^  as  !.n  tho  P^ctav  iidflsHe,. 
the  ^•'feedback  ttm  can  be  modified  to  be  -#<•'  _  /£y-r 

vt-T  ’  ^  . 


Since  jq 

V 


r 


TP‘#- 

/■ 


V/ 

‘'O.-tf— 


’«  o<. 


This  results  in  ji  coh-rol  ecraation  of  the  fosia 


ic = ■ 


c<r. 
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2*2~12 


The  closed  loop -transfer  f  unction  will  be  modified  slightly  and  be  of  the  f ova 

,  *’  •  ’  ■*  *-  or*  ^  . 

/■  ...  v  -  ;  yr" '  ■  - 


s* 


Is  VH£.  '  (. 


K,  __ 

.  VgVJT 


•  /ivr 


/^r 


With  this  form  the  static  gain  ±3  approximately  unity  at  all  flight  condition 
and  for  any  angle  of  attack. 

Itype  A  and  Type  -8  systems  are  not  practical  for  the  Tartar  missile  since  trie 
variations  in  41  idtli  angle  of  attack  is  prohibitive  a  nd  the  airframe  is  unstable 
for  large  angles  of  attack* 


Type  Caid  Type  D  are  possible  systems  for  the  Tartar  missile <>  These  are  the 
only  systems  that  will  be  considered  in  the  following  sectiohSo 


2*3  Stability  Analysis,, 

If  the  systems  were  complete  as  presented  in  the  previous  section  the  only 
analysis  that  would  be  required  would  be  an  examination  of  the  roots  of  the  characters  fci 
equation*  For  the  type  C  system  this  means  the  determination  of  the  roots  of  a 
quadratic  equation,  and  for  type  D,  the  determination  of  the  roots  of  a  cubic  equation* 

In  general  if  the  real  parts  of  the  roots  are  negative,  the  systems  are  stable* 

There  will  be  adequate  stability  if  the  ratio  of  the  real  part  to  the  imaginary 
part  is  the  tangent  of  an  angle  less  than  60s*  This  corresponds  to  a  f  y>.  *5* 

sr 

or  the  complete  system,  however,  such  an  analysis  is  inadequate*  The  affects 
of  the  following  factors  must  still  be  determined,, 

(1)  *  Instrument  responses 

(2) *  Servo  response 

(3) 0  Rate  limit  on  the  servo 

(lj)o  Expected  tolerances  on  the  instruments,  circuitry  and  gains* 

(5)o  Additional  filters  for  noise  or  bod/  vibration* 

The  most  satisfactory  method  of  examining  these  effects  has  been  to  examine  the 
open  loop  transfer  function  by  Nyquist  plots*  For  multiloop  systems  in  general 
the  open  loop  function  will  vary  with  where  the  loop  is  opened*  The  open  loop  transfer 
functions  for  the  too  systems  will  be  derived  first*  They  will  be  derived  for  the 
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following  case 

(Do  Loop  opened  at  the  wing  servo 
1*).  Loop  opened  at  the  rate  gyro. 

(3)o  Loop  opened  at  the  accelerometer. 
See  Figure  3  (a) a  >(b)  and  (c). 

For  type  C,  caso  (l) 


i  '  <$/* J  A-  4^  V  3fJ  V-  a;  4* ¥  %)] 

t  '  f 

Substituting  for  ^  and  from  equations  2„1~8  and  2v<,l-p 

'■  5,3~'  kn  ' 


/  ~~  &  5 

C  ^  *?T 


For  type  C,  Caso  (2) 


~  GJS) '*■>.(  % g/j)(% Tg/s)#,  J 

Since  negative  feedback  is  not  implied-  for  any  of  the  inner  loops  considered 
the  negative  sign  is  put  on  H  (S). 

^ _ §£L  ..  1 

f'wtfi  ■  N  ^  7  .2.  1  .  7  . * | 

~4  7  '/*£*')  *+  -&£  I 

"  wj  <r~  /  / 

4,  -  -M/*r <$*■)(„ m _ _ J 

7^/-  <2  6rs«  ^3  )  .  J 

K,  $/*)(%)/ fas,) 

!-  — 

S’  3  -  3  /•?<,  _ <r</s  j  .dar;  £  .Aisiz.. )  6^s> 


<*>> 
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Assuming  that  &A)  includes  the  integration  term,  the  open  loop  transfer  functions 
for  typo  D  will  be  the  safe  as  for  typo  C  if  K ^  is  replaced  by  Kg  <•  S0 

The  information  that  can  be  obtained  from  examination  of  these  open  loop  transfer 
functions  are  the  following.,  *  - 


For  Case  (l). 


Effect  of  servo  response  and  variatioi  thereof *  Nonlinearity  can  be 
considered  by  use  of  Johnson5?,  describing  function,' sD 

For  Case  (2)° 

Effect  of  rate  .gyro  response  and  variation  thereof*  Effect  of  additional 
filter  following  tho  gyro. 


For  Case  (3)o 


Effect  of  accelerometer  response  and  variation  tr,erc.of0  Effect  of  additional 
filter  following  the  aeeelero  eter0 

For  purposes  of  seeing  how  the  Nycuisl  plots  look  qualitatively s  the  instrument 
anc.  servo  response  can  be  made  ideal* 

(: yr/l< )  «  / 

/i)  "  1 

(5)  -  /  fa  Ti/pe  C 

-  .  Fer  ® 

0 

The  following  approximation^  car.  be  .naue  for  the  aerodynamic  coefficients  * 


<S* 

& 


ft  tz  -  SC  Q£ 


Tno  open  loop  transfer  functions  will  then  be* 
For  the  type  C  system* 


?.*>6 
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F>?r  the  type  0  .system* 


2o3“S 


203-’9 


4-  -  /9£  r  j  -t-k,  Tgifcrl 

4in  C  L  <■/  ~£  5  -.i1)  5  :  .  J 

X  -  -  (*,■«•*>)(■# )  t^£M 


•*  V"  Vm.  <  /  <9  '  P  \ 

.  .  „  c  c.  T  /^rrK>  J 
.«  u  Sj£.  )  &£ 

— . J^ai  _ ... _ _ _ 

5  -  #  *  -  #  -f-  % }  *w 


Since  negative  feedback  v.as  not  inpli nd  in  the  derivation  of  the  open  loop 
transfer  functions ?  the  characteristic  equation  is  of  the  form 

/  JS*  **  o 

This  would  mean  that  the  encirclement  of  *  1  on  the  polar  plot  snould  ba  examined* 
However,  since  it  is  more  conventional  to  examine  the  encirclement  of  •*!,  all  the  open 
loop  functions  should  be  multiplied  by  -'l  before  being  plotted*  The  characteristic 
equation  will  then  be  of  the  form  /  .j.  j  ^  l-t-Cf/s)  ^ O 

The  simplified  form  of  these  open  loop  expressions  are,  in  gens  ml ,  of  the  form* 


—  -*«■ 
%* 


*  K 


>  <?,  ;  ~M'7  j  V" 

/ ~i~  Cs  .S  ~j-  Cut  £*'r  Cy  $  ' 


and  the  Nyquist  olot  can  be  computed. ritn  little  effort..  The  casec  where  the  aero** 
dynamic  coefficient  C  is  positive,  trier*  will  bo  poles  in  the  right  nand  plane  and  it 
would  be  necessary  to  examin*.  the  countej'ciockwisir  encirclements  of  **1  aft#)- the  ncmasr 
cf  then  5  poles  Is  determined; 

These  Nyquisb^  of  the  simplified  open  loop  transfer  function  are  extremely  useful 
in  get  tin  ,•  a  qualitative  understanding  of  ohe  eye  tern,  figure*  '6  shows  a  typical  ex® tua< 

2 Ji  ■  Determination  Of  Gains* 

The  d ctual  determination  of  the  gains  will  in  general  be  an  iteration  process* 

This  is  necessary'  because  all  the  possible  factors  what  must  be  co1  sr.dersd  cannct.  bs 
•included  in  any  one  analysis  or  simulation. 

For  the  Tartar  missile  the  factors  that,  bad  to  be  considered  vere  the  following: 

(1)<  Elastic  body  coupling. 
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;r  .  .  ' 

\2)'»  Radon®  error  slope  coupling. 

Wlth  r0ii  SyStem  in  1113  Presence.of  aerodynamic  coupling, 

(W-  Control  surface  servo  limitation. 

(5) ,  Instrument  limitations. 

(6) 0  Tolerances. 

(7}o  Effect  of  noise. 

»fwT!lv>  “J3**1,8®*  °f  ®alnS  Would  bs  v:hich  8ives  the  fastest  response  for 
sf tei  all  these  factors  nave  been  considered® 

These  factors  will  be  expl®  ined  in  -detail, 

f 

2oijol  Elastic  Body  Coupling® 

- 1 .  ,.Fi;  “f3  j*  shows  a  block  diagram  of  the  type  D  autopilot  system  vdtn  the 
aadiixonal-^lasuic  boay  loop  included®  The  loop  coupled  through  the  tdpid  body 
aei  ynamxc  responses  which  had  been  the  main  loop  considered  up  to  now  can  be  • 
nS^S°?S!°  J/iS  iS  ^e“fuse.toe  elastic  body  resonance  frequencies  are  on  the 
10  rad/sec^  rad^3ec  anc*  ',he  body  aerodynamic  responses  cut  off  at  around 

in  Appendiy5!3^1,  **ttnc^on  for  *he  ®*as'"Tc  body  responses  are  derived  in  detail 

,  ,^n  P^-i-c^ar  the  transfer  functions  an  derived  for -the  case  where  the 
bfS0?,^"?3  ela3tS  bOT1.Parti*  Ottm-M*.  equation  33 
by  tfa-  first  term  m  the  series  solution  and  the  body  loaa  assumed  negligible® 

th_  J?®  If?13  'f  ^  series  solution  represent  the  various  vibration  modes  of 
the  Ooam  arivx  only  tne  mode  with  the  lowest  frequency  is  considered®  Th-'s  *s 
reasonable  since  the  second  mode  is  at  sufficiently  high  frequency  so  that" the 
filtering  required  for  the  first  mode  will  aaeouatcly  take  care  of  the  second . 
and  signer  modes®  The  aerodynamics  loads  due  to  local  angles  of  attack  were 

muht  condt*--in»elifi<>le  8f~C'C+,  cn  ,the  rr'a5e  sh&i,t3s  or  frequency  for  a  representative 
f_xght  condi  .Ion.  Tnio  moans  that  tne  solutions  for  the  Tartar  bending  modes 

computed  tor  vacuum  condi -cions  is  approximately  tie  same  as  for  moving ~air 
conditions^  c 

transfer  functions  for  the  instrument  outputs  for  an  arbitrary  force 
inpuc  at  the  control  surface  station  were  found  in  Appendix  1  to  be 
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wher® 

*  shape  of  firot  benclAng  mode  computed  with  an  a  rbiiraiy  displacement 
of  1  in  at  i  tation  X„  . 

*  displacement  at  K  °  the  station  for  the  control  surface  hinge 
line  (in). 

fCX<J  *  displacement  at  X  s  TV-.,,  the  station  at  which  the  accelerometer  is 
located  Obi)* 

(AS.)  *  local  slope  at  X  *■  X^.  the  station  at  which  the  rate  gyro  is  located 

'•H.*,  (As*e  ). 

/'A a  »  generalized  mass  defined  by  t*\  «  x\,/iA't  <  £vhera  *  r?ass  at 

^  station  ,7*£  and  oisplacement1- at  station  £  J  (lb  in  £>cctf) - 

»  frequency  of  first  benuing  mode  (rad/sec). 

ftp  “  acceleration  due  to  bpcy  bending  at  the  accelerometer  station  (g’s).' 

V'y  *  rate  of  boay  bending  at  the  rate  gyro  station  (aeg/sec). 

fy  ®  force  input  at  the  control  surface  hinge  line  (lbs). 

S  «  structural  damping  coefficient* 

The  transfer  function  for  the  elastic  body  system  v/ith  the  loop  opened  at 
the  v.dng  for  the  type  D  system  is 


2  0 1^*3 


f  H 


-=  <&>(•  j[(Si)oA(^)(^  t  ^ 


In  order  to  insure  that  there  w  ill  be  no  instability  in  this  loop 
should  be  \  *»  o3  whenwrWjfor  an  assumed  structural  damping  of  «  o02.  This 
margin  would  be  insurance  against  the  tolerances  on  S  f  u/  ,  A(y)  r  au‘toPfI°J>' 
gains,  instrument  responses,  and  components* 
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The  open  loop  transfer  function  needs  to  be  examined  only  ato*  s.nfo  since  at 
any  other  U)  the  amplitude  will  be  lesse  ‘  ...  -  - 

Substituting  or  S  into  equation  2jj»3 

M  feL;  6M<iHl(%f*v(itL,.  ™  ± 

+  <?(h)  nxj  G^u^0)Kt  7 

2  S  r'TTST.  2  J77TT  J 

The  obvious  method  for  making  A* A  j  4.  0 3  is  to  locate  the  instruments 

'  ‘a  ft  .t 


such  that 


Jry/  *0  and  *  Qo  figure  5  chows  the  computed 

first  bending  mode  shape  for  the  Tartar  missile,,  ^ue  to  packaging  requi rements, 
the  instruments  for  nhs  Tartar  missile  are  at  the  following  stations: 

Accelerometer  •  Station  \6 
Rate  Gyro  ©•  Station  80 

With  the  instj*ument3  at  those  stations  the  values  for  the  displacements 
and  slope  are 


,2 


.  «  oOl 

f /jc-iXx 

At  the  control  surface  hinge  line  '/At*)" 
The  values  for  the  remaining  factors 
31  3k5  a&M,  I  o-l-ts 
r*  =  0^39  xtr  •:*  &  ac"- 
•'  «  *02 

Substituting  these  values  * nto  equation 


( |  :::  Oi {j  w\-j )(  fr  f  ' .  e*.;  V £<..  l/ u  j,  (,/>i  f*i  /. ,  •>'-/ 


d-  ,  C  x  ivt.J  K,.  f 
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The  terms  that  will  give  the  greatest  contribution  is  the  tern  vith.K., 

Since  ^  lOK^  ^  I&CS3  ^or  ^ar^ar  missile,,  Neglecting  all  . 

tut  the  term 

IpL)  «  (j  t*s)  ( 4C1  (<<>?**  )  S  Kj  3f-r) 

In  order  to  determine  the  amount  of  attenuation  that  is  required  from 

(filter  if  required),  these  can  be  lumped  together  into 


GyfjU'tf)  5r^CiM^3  Of,. 

J«'8 

since  Hr+.(i'w)  ^  Ju>  for  type  D0 

2»Ij-b  (  ^  Gj,  (/lt/3  )  (»038U)  ^ 

Fy-  can  be  computed  from 

?T/X^  H  8 1  A  fa  (ft*  lA*.&h 

The  value  used  for  d*;-  should  be  that  computed  from  oscillatory  aerodynamics 
for  a  frequency  of  „  For  the  Tartar  missile,  however,  the  values  computed 
from  oscillatory  aeroiynard.cs  did  not  differ  appreciably  from  chat  for  the 
stationary  aeroriynamicso 

Studies  nave  shown  that  the  Kl05  condition  fives  the  largest  (  ) 

max  value  for  the  Tartar  m.issilec  This  is  assuming  that  K~  is  programed  to  vai-y 
with  total  pressure.  '* 


A  sample  calculation  will  be  made  for  th^  Hlo5  S„Lc  condition 

JX  «*  (11481)  („?&)  (l.*?)2  (1.25) 
x 

*  Itlli 

(%n  W-  A:  Of  t'jl%)  l5.%3 

The  1L  ultimately  arrived  at  for  the  Tartar  missile  for  this  condition 
vae  K,  »  ,,175  (this  is  the  value  for  that  takes  into  account  the  small  amount 
af  p  that  the  accelerometer  senses  di-.ia  ■■■,0  .Us  being  off  c0g.)., 

f£*^«iuK=  6f  ^ 

I  '  ~  •  3 

~  its 
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This  means  that  at  to  *  3h$  rad/sec  the  combination  of  gyro  response, 
servo  response  and  filter  response  must  rive  an  attenuation  of  1/10«8  or 
-  20o6  D*3„  By  serve  response  is  msant  the  response  over  and  beyond  that  of  a 
perfect  integrator 0  The  rate  gyro  is  a  second  order  system  and  the  state  of 
the  art  is  such  that  its  natural  frequency  ahe'  damping  ratio  can  be  specified 
to  any  arbitrary  values  within  reasons  *or  the  Tartar  missile  the  specified 
nominal  values  for  the  rate  gyro  are 

natural  frequency  *  2?  cps 

damping  ratio  *  0$ 

this  will  provide  an  attenuation  of -12DB  at  3h$  raa/sec.  The  additional 
-8  DB  can  be  obtained  fdther  from  the  servo  or  an  additional  first  order  filter. 
What  must  be  considered  in  the  process  of  determining  how  to  obtain  the  desired 
attenuation  is  the  ihase  lags  that  will  be  introduced  at  the  gain  cross  over 
frequencies  for  the  Nyquists  of  the  systems 

For  the  Tartar  autopilot  system  the  g;in  cross  over  frequency  for  the 
Nyquist  with  the  loop  opened  at  either  the  control  servo  output  or  at  the  rate 
gyro  is  about  liO  rad/sec „  This  can  be  seen  by  plotting  the  simplified  open 
loop  transfer  functions  or  by  making  the  assumption  that  at  the  frequencies 
concerned 

Jx.  Pd  g 

^IA  A* 

and  therefore  Jue  gain  crossover  is  at  Lo  ^  EK^o  The  simplified  Nyquist  for 
*»/, .  Ml.S  S.L.  is  shown  in  Figure  6,  * 

At  JUO  rad/sec  the  additional  phase  lag'  introduced  by  the  ~?te  gyro  nnd 
the  required  first  order  filter  are 


(1). 

From  the  rate  gyro 

Hi5 

(2). 

From  the  filter 

3Jf> 

(corner  at  2li  cps) 

The  system  must  be  able  to  tolerate  this  additional  phase  lag  a  nd  still 
have  an  adequate  phase  margin  of  cz  30°. 

Z£  it  turns  out  that  this  is  not  the  c ass,  the  possible  fixes  are: 

(1) o  Move  the  instruments. 

(2) .  Design  a  complex  filter  that  gives  the  req  dred  attenuation  at 

but  less  lag  at  system  gain  crossover-  frequency. 

(3) .  Find  tome  means  to  redice  K^o 
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For  the  Tartar  missile  the  additional  28*  of  phase  lag  at  IV  c  bO  rad/sec 
was  tolerable,. 

There  exists  therefore  a  choice  of  putting  the  first  order  filter  after  the 
instruments  or  incorporating  it  in  1h^  transfer  function  of  the  control  servo® 

If  the  filter  is  placed  at  the  output  of  the  instruments,  the  servo  must  be 
designed  to  introduce  as  small  an  amount  of  phase  lag  as  possible  over  toe  nominal 
PC'5  at  IjO  rad/sec  o  Since  this  woulo  add  to  the  28f  already  introduced#  The 
additional  amount  that  Tartar  can  tolerate  is  on  the  order  of  ®  The  disadvantage 
of  this  is  that  the  requirement  on  the  servo  is  prohibitively  stringent®  The 
advantage  is  that  if  the  requirement  can  be  met,  it  is  possible  to  use  high  gains 
in  the  roll  system  -without  complex  networks®  This  comes  about  because  the  roll 
system  uses  the  same  control  surfaces  as  the  autopilot  and  phase  lags  in  the  servo 
are  introduced  in  the  roll  system  also® 

The  advantage  in  incorporating  the  filter  in  the  servo  res uonse  is  that 
the  requirements  on  the  servo  are  deliberately  relaxed® 

Both  methods  are  being  considered  for  the  Tartar  missile® 

2oh<>2  Radome  Error  Slope  Coupling® 

The  equations  which  show’  toe  existence  of  the  coupling  due  to  radome  erx*or 
3lope  will  be  derived  first®  The  angles  squired  for  the  analysis  are  shown 
in  Figure  7® 

These  angles  are: 

<f  <*  Angle  between  line  of  sight  to  target  and  reference. 

0  »  Angle  between  seeker  centerline  and  reference, 

Y  *  Angle  betv/een  missile  centerline  and  reference® 


/3  *  Angle  between  see  ker  centerline  aria  missile  centerline  - 
t  «  Angle  between. line  of  sight  to  target  ana  centerline  of  the  seeker¬ 
s',*  mirror  angle  due  to  rc.doms  refraction 
j3  ft  ■  Look  angle  °  (  <r  ~  Y  ), 

The  relationship  between  these  orgies  ai'e: 

£  -  <s~  Q  —  y-  V'-/B 
®  ~  -h  f 
~  £  +•  £f 

£  will  J.n  general  be  a  f  met:  on  of  the  look  angle  a  +  £. 
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£  t:  {■(fii-lj 


and 


if 


<*£i  "  j. £j L  4 CjSt+fz} 

dC^^-e;  _;<■- 

jj is  assumed  constant  and"sei  equal  to  jt* 

C/9+i.) 


2oltvo  £r  ~si.{  j&f-Bl small  regions  of  look'  angle  this  approximation  can  be  rnadG,. 

Since  the  dynamics  of  the  seeker  is  not  required  for  the  analysis  in  this 
section,  it  will  not  be  presented*  It  is  sufficient  to  state  that  the  function 
of  the  seeker  ie  to  continuously  track  the  target  and  in  the  process  of  doing 
this  provide  informtion  re,  arding  the  rate  of  rotation  of  the  line  of  sight 
angle  to  the  target  (  $his  information  is  immediately  available  since 

if  the  seeker  tracks  the  target  with  a  reasona'  le  Degree  of  accuracy,  the  rate 
of  rotation  of  its  centerline  in  space  is  identically  the  rate  of  rotation  of  the 
line  of  sight. 

The  ho’dng  guidance  equation  i3  ideally 
~  <r~ 


vx 


where  -A-  is  a  constant  ^  y  and  V  is  the  relative  velocity  between  target 
and  missile,  Gg  (s)  is  the  transfer  function  of  the  guidance  filter* 

for  the  Tartar  missile  <r~  uiroi'. action  is  obtained  from  i  or  the  head 
rate  'gyro  signal  plus  the  rate  of  the  tracking  error  signal* 

From  equations  2»l4-£  and  2<h~6 

&  1-  i  ' ~  6  E  +/!,(  £  +  &) 

i  * 

•z  (  i-t-svj  a” -si..  ^ 


and  since 


h  c  (  3  W  ~  V/»  r  { 


l?VSr 


n  C  ~  [  (  i  h.'i  )  <r-  -  A-  •>  ]  Cr^  ( A .J 

j  y  >fj~  * 


For  the  Tartar  missile 
Og  (s)  *■  1 

t r*iz&rvn^®s) 

This  was  deterrined  by  &  study  on  the  amount  of  filtering  rsqud  *%d  /or 
adequate  homing  guidance. 
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t*  £-A„/  ,--  ■>"  *  -  ^  'Vi'  *?  ,  -  - 
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The  complete  equal-ion  for  is'  therefore 


This  shows  that  there  is  an  additional  ^  feedback  into  the  autopilot*  A 
block  diagram  of  this  additional  coupling, loop  is. shown  in  figure  |U 

The  open  loop  transfer  function  with  tire  loop  opened  at  i4*  is 


n/ric 


j  £  <4&* 

The  simplified  open  loop  transfer  function  for  n/nc  is  of  the  £0111; 

i 


a  I;  tk$  -u 

^/n  is  of  the  form 


n.  -  t'.****  c  i  -h  %  .w 
vl  r 

other  e  T.tp*’ 


(  frot+y  sccktt'A  2. 1  J 


5 


/?£-£?  C 

The  open  loop  transfer  function  becomes 


V'i 


/I-  _ 

/ 1  f ,  ^6^  '){(+*  2uv,dJb  ,*•  §  '&*•* 


When  is  .jooitire  there  is  negative  feedback  cn  the  loop  and  when  it  :in 
negative,  positive  feedback*  Since  1  ith  positive  feedback  this  .loop  becomes 
extremely  difficult  to  sfcablilisa,  fU  can  be  biased  so  that  ir  effect  it  .ms 
only  positive  values-  This  is  cone  0/  adding-fo  a  Vy/ some  paid  rate  gyro  signal 
*  Jl>  urn:  y>  f  * 


then 


4L  : 


.  s*  y&SJjL  J.l%k-dzL 


'«  V“M  £  /  ■**.  r.k)tf  *■ ,  k  trdj 


The  value  of  ( ;v  /v  )  will  never  be  negative* 


tXTStiT 


* 


( 


<  Xv.\  \ 


pnr 
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V;v  />.  ?  v:F.pr  raS  #•/»*£  “  *'  oOh,  the  quantity  (  )v£  <A-  ^  o£ 

'  JL  £*• 

trhila  ’"7”£.  J;  which  equals  ffgZgC;  *'■ can  be  ^  3arRe  as  &*°  a't  high  altitude 

'■  conditioua.  where  A  decreases®  The  autopilot  time  constants  TA  and  Ta  must  be 
adjusted  to  stabilize  the  loop  for  these  conditions®  At  low  altitude  conditions., 
there  is  no  stability  problem  in  this  loop  since  is  small  and  the  attenuation 
from  the  guidance-  filters  is  more  than  sufficient®  Figure  9  shows  a  Bode  pilot  of 
the  tw/4'r&,  open  loop  transfer  function  for  the  M  £<>0  h  **  50*000  ft®  flight 
condition.  For  a  fixed  flight  condition  varies  with  «£  but  T/(  varies  roughly 
^||te>.8afe_  iwni^  '%ue  ijgaAgtajLw^pg:  the  same  degree  of  stability. 

An  alternate  method  for  stabilizing  this  loop  .is  the  following;  Design  the 
auto  pilot  for  maximum  possible  speed  and  design  the  guidance  filter  to  stabilize 
the  loops  The  advantages;  and  disadvantages  of  both  methods  are  discussed  in  the 
section  on  guidance  computer. 

Sol Compatibility  With  Roll  System  In  The  Presence  Of  Aertxiynsrile  coupling 
Derivatives,, 

This  subject  :1s  covered  in  full  in  the  section  on  roll  system.  It  is 
sufficient  to  say  at  this  point  that  the  gains  must  not  be  determined  without 
consideration  of  the  roll  yaw  coupling  problem.  In  general  the  roll  system 
alone  cannot  be  designed  to  stabilize  extrema  roll-yaw  coupling  instabilities. 

.'light  modifications  of  the  autopilot  gains  can  improve  the  situation  considerably 0 
As  an  example,  low  damping  ratios  /  in  the  pair  of  complex  roots  of  the  autopilot 
is  detrimental  in  the  roll  yaw  coupling  loop<>  This  low  $  can  result  from  the 
fact  that  the  stability  derivatives  for  yaw  motion  :is  different  from  that  for 
pitch  alone  when  the  missile  is  pitched  at  a  large  angle  of  attack. 

2ohoU  Control  Surface  Seivo  Limitations. 

The  following  limitations  ’"•ust  be  considered. 

(1)  Maximum  practical  servo  response® 

(2)  Maximum  wing' rate. 

(3)  Effect  of  wing  loads  (.hinge  moment}  on  (l)  and  (2), 

*  {!' )  Non  linear  phase  lags t. 

It  has  already  been  pointed  out  in  section  2Jui  that  the  consideration  of 

*  the  maximum  practical,  servo  response  determines  the  method  by  which  body  bending 
mode  coupling  loop  is  stabilized® 

The  effect  of  maximum  wing  rate  on  the  stability  of  the  system  can  be  detemuned 
by  use  of  the  describing  function  on  the  Nyquist  of  the  open  loop  transfer  function 
for  the  loop  opened  at  the  swing  servo  (  )„  This  will  immediately  snow 
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that  the  type  C  system  with  the  position  feedback  servo  is  more  susceptible  to 
non  linear  oscillation  .than  the  type  D  system  vath  the  rate  sax'  p.  A  oossible 
means  of  improving  the  situation  is  to  add  another  non  linearity,  such  as  a 
frequency  sensitive  limit  on  the  input  to  tho  servo* 


The  effect  of  wing  loads  is  to  induce  the  servo  speed  of  response  ana  possibly 
reduce  the  raaximu '  wing  rate*  Tnis  must  be  considered  as  an  extreme  tolerance 
on  the  servo  response,. 

The  non  linear  phase  lags  iromtceau  space,  sticky  value,  etc*,  should  be 
estimated  and  considered  in  the  (  ht/t  )  Nyq  '1st  plot* 

The  control  surface  servo  transfer  functions  are  derivec.  in  section  3e.5»2 
T  is  is  the  section  on  scr.o  limitations  that  affect  the  roll  system*  The 
derivation  is  in  this  section  since  the  servo  response  is  more  critical  for  the 
roll  system,, 

2oUo5  Instrument  Limitations. 


The  following  factors  must  be  considered,, 

(1)  Compatibilitory  of  specified  dynamic  response  with  the  dynamic  range 
required. 

(2)  Compatibility  of  the  dynamic  range  required  with  the  null,  noise, 
resolution,  g  sensitivity  and  linearity  characteristics 0 

The  desired  dynamic  response  can  be  determined  by  the  method  suggested  in 
section  2o2j<>1  if  the  response  is  to  be  used  to  stabilise  the  body  bending  mode 
coupling  loop*  If  it  is  physically  possible,  this  is  desirable  since  the  phase 
lag  from  the  instrument  need  not  be  designed  out  of  the  system* 

The  effect  of  null  unbalance,  noise,  and  g  sensitivity,  can  be  computed 
by  assuming  that  these  are  the  only  inputs  to  tue  system  and  computing  what  the 
output  (N?t)  will  be  with  these  inputs* 

As  an  example  the  effect  of  additional  inputs  at  the  rate  gyro  on  steady 
state  condition',  fill  be  computed* 


The  control  equation  is 

Ar-l  v  ft,  (  A  ~  A$J  ,&,)  (  f  f  iSJtf  «*) 

where  ceuj  be  the  null  value,  or  the  amount  due  to  g  sensitivity  ;  -  'W  ,/•>) 

For  steady  state  conditons  ^  4- 
K  ;  (  -n  -  t\  £  ]  ~ 

i)i%'  K  HK' 

M  *,3 


O  c 
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This  gives  the  amount  of  error  introduced  by  the  rate  gyro  in  the  steady 
state  condition* 

2oho6  Tolerance* 


The  gains  nr.st  be  chosen  to  be  able  to  stand  whatever  tolerance  cun  be 
practically  met*  For  t  he  Tartar 9  the  gains  were  chosen  to  give  adequate 
performance  with  20p  variation  in  any  combination,  'lids  was  based  on  the 
expected  accuracy  of  the  P,  meas  ring  and  gain  setting  device  and  the  component 
tOlfi  WttCfifit  T, 


Effect  Of  Noise, 


The  effect  of  noise  from  any  possible  source  can  be  computed  by  assuming 
one  source  at  a  time  as  the  only  input0  The  cutout  considered  diould  be  both 
(  n.  )  and  i  o 
r 

As  an  example*  ji.f  there  is  noise  at  the  acc.lero.ieier  and  it  is  desired  to 
find  its  effect  or.  3  the  transfer  function  between  2',  and  must  be 
derived*  -^he  control  equation  with  a  noise  input  at  the  accelerometer  is 


Gpl*)  if,£n  f  nJ 7 


where  (S)  is  the  filter  following  the  accelerometer. 


Since 


r’  •{  ^ )  &.Js>  4  (■>!  i-c, 

'f‘  ~(^-J  4 (S)  CrJ^  ic 

(  A>)  Cpr2Sj#+'~&) 


Ct.%  *  r3xt 4  J  ic 

^  •*—— . . . . . . . . .  . 

'Wo*  1  ~  %.(*>(&)  -  (%  t 

%Gu) 

can  be  computed  iVoro 


The  RMS 


fi 


A, 


Hi 

air* 


'S» 

(jr  { M')  $  ' 


X) 


(*s 
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where  is  the  spectral  density  of  the  accelerometer  noise  in  g^/rad/sec. 


of 


Normally  it  is  sufficient  to  throw  out  all  the  terms  in  -the  denominator 
except  the  1  and  G  (t>)  becomes 

£Pu->  *. 

Tills  is  justified  for  two  reasons  (l)  the  feedback  terms  are  attenuated 
quite  heavily  and  (2)  the  RMS  calculation  is  more  conservative  without  them* 

» 

The  can  now  be  simply  comouted0 

3o  sr.r.xx  dpphst- 

3ol  Description  cf  System 

Appendix  II  describes  the  system  that  is  to  be  analyzed,.  This  system  has  the 
following  degrees  of  freedom: 

(1)  y  a  translation  in  the  yaw  plane* 

(2)  W  *  rotation  about  the  c.go  in  the  yaw  plane. 

(3)  0  *  rotation  about  the  missile  centerline* 

(Jj)  ly  <*  control  surface  deflection  to  produce  yaw  f oi-ce  and  moment. 

()>)  <>  ®  control  s  irface  deflection  to  produce  rolling  mo-  ent. 

The  angles  said  the  stability  derivatives  required  are  all  defined  xn  Appendix  II* 
30g  Form  of  Control  Equation. 

The  type  D  control  system  will  be  assumed  for  tea  autopilot. 

The  roll  system  control  equation  can  again  have  any  number  of  iorms#  The  form 
that  the  Tartar  missile  uses  is  the  following 

j.M  f.  -  -  bte£c£x£*ii2J 

When  the  equation  is  actually  implemented,  it  takes  the  form 

3  2-1 '  Sc. -  k,(4>'+k>  4>')0  /-  rs) 

Figure  10  shows  a  block  diagram  of  the  system.  0  ^  (S)  is  the  fz-se  gyro  response, 
G  ‘j.  (S),  the  rate  gyro  response  and  p  and.  ^ ' the  outputs  of  the  two  instruments, 

3*3  Simplified  Transfer  Functions 

The  input  to  the  roll  system  can  be  either  a  command  to  roll  or  an  extraneous 
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•#  , 

roll  torque  (  ^  ) .  The  output  or  the  controlled  variable  Is  0.  The  simplified 
closed  loop  transfer  functions  with, these  two  inputs  can  be  derived  for  the  system 
by  assuming  that  the  roll-yaw  coupling  derivatives  are  all  aero,  the  instruments 
have  perfect  responses,  and  the  servo  is  a  perfect  integrator^ 

The  closed  loop  transfer  function  with  a  .c  lied  for  roll  angle  input  will  be 
derived  first,. 

Using  the  form  of  the  control  equation  3»2~1/ 

$  »  x*(i+rs\(4t--4)-Kt,i(to+rs)0 

Where  (d-<£)ia  the  free  gyro  output, 

%  -  [K4  +  (KT  +  fG<s') S  +  KTS* ]<j>  +•  K U*TS)4>t  ■ 

The  roll  moment  equation  is 

f  $  +  G  6 


b  (9 
g  "  s(s-F) 

Substituting  into  the  control  equation 


'  s'Vs-r) 


V-  A;  *  t  +  k0  r4 )  s  ♦  Airs  7 


*K<(i*Ts)l 


and  therefore 


3«3~1 


/  7  5 

i  +  0(.T  -£  )s%  +  (r  +  Ka)>  +  l 


For  a  step  (p  the  roll  morent  equation  is 

$  *  6  S  r  F  <p  +■  fa 


l  * 

~  6  c)  -  0c 
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with  the  control  equation7  -  * 

!.  i  +  (K*T  +  K*  K„h  +  tu  K4  75*14  *  s  £  *  ° 
* 


fT  ’ 

rO. 

and  therefore 

4> 


3„3"2 


5  3-  SZF  *  6K4  6(c4  T  %  )&■-+  &  tf  n  r S* 


.s  . 
<5  A4 


+  (ti,7- Ifa)  5X+  (r+-ic,)5  +  i 


Equations  3o3“i  fcTjd  3°3-*2  give  the  basic  characteristics  of  the  system,  The 
roll  system  in  the  absence  of  voll«yaw  coupling  is  «  third  order  system.  The  static 
gain  for  a  step  S0  is  iaentically  unit y0  With  a  step  <p^r  the  system  benaves  like  an 
imperfect  differentiater. 


There  are  again  tnree  possible  places  to  open  the  loop. 

The  open  loop  transfer  fu  ction  for  the  loop  opened  at  the  servo  is 


3«3-3 


L  (Kt  +  feS  *  A*6  **06 


<>  in 


S\:S 


r) 


3 oil  Analysis  of  The  Complete  System, 

The. equations  of  motion  for  the  five  degrees  of  .freedom  are  {from  Appendix  II) 
A/t  f  Ely  4-  -  JS  ~  0  ,&sVl ,  Q(  $  •/'  Cd-iLt  « e3 

C  6  v-  £  U  4'  M 


when? 

rearranging 


iy  -  k; ’(A/i  *  *>/  +  M 1  >  *  K>-  y  +  •  *5  y' 

<$i  +  F  i  +  tij&  +  l  /  - 

Sc  *  -  {7<*  ?*“  Ks  s  r-  ’l;\  .5 7'^)  <p 

K1-  JL  u 

'  184s 

{'/?-$)£  +  3  ly  -  Co-xl  •fti-o  y-  V-  <^'1* 

*  « 


c  /» 


V  c  ?I 


a,  y  -//tff  ~  0 
V-  V  «f  -  0 
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(equation  corit- nvid)  /) yS  f  * g  -  5  )ly  +  (KZ*K^  X  *f  +  M K%  T 

~H£  -  £)<§>  —  (•?  ^ 3?  o 

The  term  ■yt  is  omitted  from  these  equations.  Since  for  the. :T^rtar 

it  v&s  found  to- have  negligible  effect.. 

these  equations  a  re  of  tne  form 


where 


4/ 

+  Art  ly  f  P>  IS  V*  ~ 

- „<u^  «.  ^  -  Af  f 

fa  ft 

^  *lt  ?,  *  fa  y-  * 

-htt 

fa/1 

*  ^31  * 

-  />-!  K,  0 

ftp 

-i-  L  ly  - 

0  *  ft -S 

ri\ 


$  £4  (P  +  fi$.£  c 


$rs  -  -CtjcL  «o 
■«»,  =  c 

HfS  :!  ~5 

4,,  - 

-  k,‘b-s 

fa  =  *1  <•  *5  * 

rf*.  --  s-f 

,flds  s  -<s 

'  Kq.  t  /vs*  S  V-  S' 

Afs"  =  s’" 
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3oka2 


defining  the  determinants 

— r-  :4irr  r-.- 

^Zl  Az.z  $23  * 

^13  : 

A34  ^ 

A**  Aj-5* 

ar>4  2- s  Ha  t  flu 

3oli«l  ' 

$  »  1 

.,  .  .  ^W 


S  *» a>-4 

4^” 


3oli**3 


^  < 


At  Aid 


(P  ~ 


4t  ^3 

^z.t  Aj 


3oll-il 


ijp  * 

/ 


{Vi  $i»  #»3 

4.  ^5  *  ,  /l/  fl« 


Tl 


r  u 


•h 

Substituting  aquations  and  3^*2  ~~ 


»  !h 


Ur'  x%  ^ 

m’  Ax  fhi 


atjt,  A 


At 


aJ 


3oli-»5  *7 


T 

V 


£ 


'3 


■V  /? 


f  ^ 

:  Ha 


^3 


/? 


35 


f  L 


fi*$  \( 

1A>  $  gg:  0co 


A 


'3 


A 


function  reTtaTI^ute'sSfl'.J  "®:'d  sld"  °r  ein-it-.cn  J.1..-5  )..-.  t„, 
«*-« 11  t,w  1,1  ■ -* 


opc?n  loop  transfer 
0  Me  rc).j  s  at C-r., 


bqiation  3oi»«S#  can  be  .m  in  the  form 


{  -  H  iv,,  o„  h>t 

j 

✓ 

1 

1 

1 

1 

^J3 

l,  <Iv It  JSV, 
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~0 

vi/  4 


+RAV/~ 


Substituting  the  expressions  for  ,.4-*  and  expanding  the  determinants 

•  *-) 

_  sUs*  Of  O  S?*~  <f<F Kj  f  )s  —  £f  K?_  *JJA 


^  -  /4n«v U  Cio  ^  -Aw»*  A*!  $•  )s  -f~  .CA/trtS  &io  C  ffjd 

<1SXV 4  -  -  ,-  I  II  r  -  «■  -*|  ■  I  I  I  1-1—  I  II  I.  n  .ii-t.i  I-  -I  i  -  ii  ~*  -  -  I  -  T  - ■-  - 1  r 


dim*' 5. 


jdt-  ••  $ -s 7  ftl  £ 'K-j  'h  A/  c<i^  (ko  ~  f'/B  iXi)s  -  ($ /V 8 )(*' , !CiC cx o  -/  ) 

^  *  ’  &\s 

~  '  .  77  ~  “ 

3  :s  (A  -y  B Kt'  -b  f-  K%)  5^  y- [c  "/rfe.  *  Aj (flP-Bi '-)]$ 

-y-  (%,  Cvo  oi.a  y1-' Ay,) 


(/> 

)Ap>  +  L  ?y 


S1*  •+•  f  Cy  .fjT/  -» )  S £7  /fg-  S  -f* 


S  **•  /(V 

? r^lo7*~6KFr~ 


+  C?  Kd 


''-■\2ar  H 
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If  tee  encirlemeiit  of  *i  on  the.  Nyquist  plot  is  to  be  examined.  ~  :*jfy  is 
plotted  since  negative  ..feedback  is  not  implied  at  the  point  where  the  loop  is 
opened,  '  ~  -  • 

Of  the  two  terms  on  the  right  hand  side  of  cq nation  3ou~6  the  first  is  referred 
to  as  Ihe  0  loop  transfer  function  and  the  second  the  £  loop  transfer  f unction,, 
Examination  of  the  Nyquiit  plots  for  these  two  loops  and  the  vector  sum  of  the  two 
give  an  indication  of  the  gain  margin  exist, in;  for  the  system  for  variations  in  H  and 
L  top ether*  Ho  information  regarding  the  phase  n argin  in  the  %  loop  alone  or  the 
<Ly  JLoop  alone  can  be  obtained,,  In  order  to  obtain  these  information  the  expressions 
for  bo„.  and  tfe*.  must  be  derived,,  This  can  be  done  in  a  manner  similar  to  how 

$v<x  v>  - J* 


the-expression  for 


was  derived. 


Figure  12  shows;  the"  Nyq  list  Hots  for  the  <f>  loop,  the  £  loop  and  the  sum 
of  the  two  for  the  Tartar  missile*  The  flight  conditions  are  Pi  £  ,  30,000  ft. 

altitude  and  *  ?Q°„  This  was  found  to  be  toe  most  critical  set  of  conditions  for 
the  Tarter  missile..  The  roll  system  gains' t  ere  chosen  to  attenuate  the  amplitude  in 
the  region  of  It  to  ItO  rad/sec.  For  more  extreme  conditions  (large  .©<  5  oigher  Kach 
number/,  the  following  happens:  The  zero  frequency  amplitude  of  the  $  loop 
increases  and  approaches  «1„  The  mag  itude  of  Ktj  anu  Iy  required  to  attenuate  the 
amplitude  sufficiently  in  the  It  »  itO  rad/sec  region  becomes  impractical,,  The  zero 
frequency  gain  of  the  £  loop  depends  only  on  the  aerodynamic  coefficients  and  so 
if  it  approaches  -1,  there  is  no  sim.-le  fix. 


3.5  Determination  of  Gains 


(t) 


or  ;  are  not  critical  since  the  guidance  a  nf  ornate  on  does  not-reqt’irc  Ine  roll 
attitude  of  the  missile  for  a  reference*  This  means  that  the  main  criteria  for  the 
system  design  is  adequate  stability*  In  designing  the  system  for  adequate  stability 
the  following  factors  must  be  considered, 

(1)  Tioay  torsional  mode  coupling 

(2)  Control  surface  servo  limitations 

(3)  Instrument  Limitations 

(It)  Tolerance 

(50  Noise 

In  general  it  has  been  found  that  Lne  method,  for  stabiliai <g  the  system  is  to 
increase  tee  gains. 

The  manner  in  which  these  factors  affect  the  maximum  gains  usable  will  ’c 
explained. 


£  ,1  pm 
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Jo^ol  Bcdy  Torsional  Mode  Coupling 

The  same  equations  that  ler*  used  in  the  bending  mode  coupling  analysis  are 
applicable  here  except  that  the  torsional  mode  shape  is  used  and  the  input  is  a 
generalized  moment  input  at  the_.tailo  - 

The  rata  gyro  output  can  be  computed  from 

<f>  1  «  JOIJl. _ , _ _  ’  £1^1  /V)  ^  ' 

.  ->  IS  . h  .  Ti  '  b  -j.  Y  /  JT 


ST  +  2$uJTZ-h«i'r  ’'  ' 

where  J  «  generalized  moment  of  inertia 

compirfoi  angular  displacement  at  the  rate  gj  *o  slation. 
computed  angular  displacement  at  the  tall 
'  //If  »  torque  input  at  the  tail 
lUf  a  torsional  natural  frequency  (first  mode) 

The  same  information  can  be  obtained  from  vibration  test  ciatac  The  test 
data  gives  the  following  information  to  the  STV--6  missile,,* 

9?S  in  lb«  at  the  tail  produces  a  maximum  of  151} « 2  -x.d/sei/  angular 
acceleration  at  all  stations  from  80w  forward  i  hen  the  frequency  is  Hi?„5  cps0 

This  say3  ,  " 


m  f  -  (% %)K  (» &) 


(V„)  <t> 


<ysf}2  /rfT  (-ft  Hi) 


/2/3 


Hitt  • 


The  maximum  open  loop  gain  for  the  torsional  mode  coupling  loop  becomes 

£,,  .  /I  i/JlV  Mt\/  /.s 
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since  K(,  ^  K?  ($&&>) 

■A}  ~  (fcf  -f-K^)Qs(s)  *  ,^4 

nv\A  f 


I 


$■ 


w, 

/) 

<«M 

/«• 

/ 


-r  r  l4Si  )\  S  d  ft? 

£ 01*  ft)  fS  >  S.L  f  M  ~ 

S7/  /’t'' 

/t)  --1  (zr(ttur)('6ii)(-n<$(s<1')  3  -?4<V'wa} 

'  r“”  •**$#.  /V 

This  says  that  at  this  condition  the  attenuation  require  J  from  G,,(da/f/.or 

value  of  *  3  for  (  )  wax  .is  '  v 

6fQwr>  ~  rl  =  '?7i 

-  •  •  - OB  ’  "' 


, .? 


*1  *  ■  o  (5 


A  rate  gyro  with  a  natural  frequency  of  85  cos  s  nd  has  been  specified 

for  the  Tartar  missile*  This  will  provide  an  attenuation  of  "?o5  DBo  The 
additional  «1  DB  can  no  obtained  easily  frem  the  fills*  teat  is  normally  required 
after  the  instrument  demodulators. 

The  torsional  mode  coupling  loop  therefore  does  not  present  too  great  a 
problem  for  the  K>,  that  hza  been  picked, 

3*5=2’  f  ritr;  1  Surface  Servo  Limitatione 

For  i  he  roll  system  it  has  been  found  that  torsional  mode  coupling  i;.  not 
too  great  a  problem*  It  is  possible  therefore  to  have  a  higher  crossover  frequency 
for  the  loop  opened  at  the  servo.  From  equation  3 > 3-*3 


The  cross  over  frequency  can  be  roughly  estimated  by 
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ana  u J  crossover  ^  GK^ 

For  the  Tartar  missile  this  was  set  to  be  3%  90  rad/sec®  At  this 
frequency  the  phase  lags  introduced:  by  the  servo  becomes  extremely  significant® 
It  has  been  mentioned  in  the  section  on  autopilot  design  that  there  is  a 
limitation  on  the  maximum  speed  of  response  obtainable® 

The  type  of  control  surface  servo  used  will  be  analyzed  briefly®  Figure  12 
shows  a  block  diagram  of  the  servo®  The  servo  amplifier  is  not  an  operational 
amplifier  and  the  input  network  must  be  included  in  the  analysis  of.  the  system® 

The  open  loop  and  closed  loop  transfer  functions  of  the  system  will  oe 
derived  using  the  following  notations  and  assumptions s 

Foreward  loop  gain  «  K  *  Kg  Ha 

a  771  i  % 

-  / 

Saur'i «>  i  r?\  Jdncc  -/-or  »  & 

Fee  d  <3he.  l<  po~t  ]yv,pe<ia/ic«  ~  ° 


The  open  loop  transfer  function  with  the  loop  opened  at  G>  the  input  to 
the  servo  amplifier  will  l>e  d  rived  first 
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The  closed  loop  tran  fer  function  is  computed  by  summing  the  currents  to 
the  grid  (G) . 

-f  «i.  se,  *-  e,  S'Ce  =  «,  *-SCc,*Ct)) 


f/  'i  2  5  7 ^5  y 


substituting  for  eg  in  the  first  equ  ion 


Hi  +  ±  (?0  S  «  .Zji 


t+ur^b  tliL} [# *'s (c> +c*> ] 


e,  +  e*.  s  -  -  e0  jjs  c2  2  *■  sfr  *  2  Sly  5  +T/ Sz)f ~  t 

*•  -  C’4  *+  ^  +  (  1  *)s  ^ 


-7-  ?  ■»  **" 

+*  If  ?5  f  lv/ 
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,  Tv"  (ct+c,)&  %i 


3o!>%2«2 


4+ *t* 


t  +  2 J  7p  %  s  * z T  r*  to  *4 ) »  *  T\-  2  x  Vto*  T  ^3 

x(jL*ct*)  ~'  ‘  *rt£.+cte)  ’  *  «7r*c;A 
R  K  X 


With  the  -actual  system  the  inputs, ,  is  made  up  of  K(jr>  "/”■&)  ,  j{^  '*£'  and 
j<4(4  +  4  4)  while  input  is  made  up  of  y>  and  S/fcS) 

An  examination  of  equation  3*£»2-l  shows  that  the  gain  cannot  be  increased 
arbitrarily  without  increasing  R(C(  V’  C?J  at  the  same  time  resulting  in  no 
improvement  of  the  response?  In  effect  tns  valve  natural  frequency  sets  an 
upper  .'Limit  on  the  speed  of  response  that  can  be  attained  without  compensation? 
TheorfeiicaHy  the  response  can  be  improved  by  oe  of  a  ieac-  network? 
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Poi*  the  Tartar  servo  a  lead  network  is  not,  used  because  of  hardware . 
considerations,  A  passive  lead  network  would  attenuate  the  D«C„  gain  requiring 
that  the  gain  of  the  amplifier  be  incr,  ased,  This  was  not  practical. 

The  parameters  for  the  Tartar  servo  as  they  now  exist  are 
/v  =  ^£>o<o 
C,  r  -  *f~~' 

/?  ~  • /3 

C.  2,  =  * 

With  these  parameters,  the  phase  lags  from  the  servo  over  and  beyond  the 
ideal  90  *  are  at  bO  rad/sec  U’ 

?0  rad/sec  1$  °  for  nominal  conuit.lons*  These  were  experimental ly 
determined 0  With,  loads  on  the  surfaces  they  are  expected  to  increase  by  about  90',iv 

The  Wyquist  ..'or  the  simplified  open  loop  transfer  function  (  £>  /*',  ^  )  for 

the  Tartar  roll  system  shows  that  there  is  a  phase  r.iargin  of  about  6f>°  ,  If 
the  servo  subtracts  Ssr  °  and  the  rate  gyro  response,  -as  10"  t  there  remains 
Uo  ,  The  .filter  follow!  'g  the  demodulator  which  is  required  to  remove 
the  800  cps  ripple  can  subtract  another  10*  and  the  remaining  pnase  margin  will 
be  30  °0 

If,  however,  the  other  effects  mentioned  in  the  section  2 „i>«3  contribute  more 
phase  lags  it  would  be  necessary  tc  use  a  complex  lead -lag  network  .following 
the  demodulator. 


The  affect  of  an  unbalance  ir.  the  servo  amplifier  can  be  computed  by  assuming 
ah  input  exists  at  the  ./id 


P 


where  2  lo  is  toe  equivalent  input  at  the  grid  of  the  servo  amplifier. 


then 
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or  for  steady  state  condition 
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The  feedback  ■  ain  thej'ef ore  determines  the  effect  of  the  unbalance  The 
effect  of  this  equivalent  <fbias  can  be  then  computed  for  the  overall  system 
by  substituting  it  into  the  control  equation 0  -  — -  - 

3 <5 ° 3  Instrument  Limitations,  Tolerance,  and  Noise 

The  effect  of  these  factors  can  he  handled  in  a  manner , similar  to  that  for 
the  autopilot, 

ito  GUIDANCE  SYSTEM: 

The  guidance  system  will  be  defined  as  the  system  with  the  following  input  and 
output,  .  — 

Input  «*  geometrical  line  of  sight  rate 

Out  nit  “  missile  acceleration 

With  this  definition  the  receiver,  the  seeker  dynamics,  the  autopilot  and  the 
guidance  computer  is  includede  The  receiver  mil  be  considered  only  as  an  imperfect 
error  sensing  device  and  a  3ource  of  noise.  Figure  13  shows  a  simplified  block  diagram 
of  the  complete  system.  The  angle  s  and  terms  used  have  already  been  defined  in  figure 
7  and  section  2-.li-?,  The  seeker  control  system  and  the  guidance  computer  will  be 
discussed  briefly, 

llol  Seeker  Control  System 

The  seeker  control  system  has  two  £:  ‘ons  (l)  to  track  the  target  i9e,  to  lafiep 

the  dish  nutation  axis  on  the  targetj  anu  to  maintain  the  wish  nutation  axis  in 
space irom  rotating  -when  their  is  no  signal  from  the  receiver.  The  control  system 
shown  in  the  bloslc  diagram  is  that  for  a  BPN~2ij  type  system. 

The  operation  of  the  s/stem  will  be  explained  qualitatively  for  the  following 
cases:  (1)  when  the  line  of  sight  is  rotating  at  a  constant  rate  and  (2)  when  the 
missile  is  pitching  at  a  constant  rate.  Figure  13  must  be  referred  to  in  order  fco 
follcy  the  explanation. 

Case  (l)  when  the  line  of  sight  is  rotating  at  a  constant  rate  (  &r  ),  the 
head  must  rotate  at  the  seme  rate  to  contina  tracking..  This  means  tnat  there  will  be 


and  for  the  Tartar  Servo 
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a  /3  equal  to  and  since  (ft  is  assumed  aero  &  will  also  equal  .  cr  „  If 
/3  is  to  be  constant  howeverthe  input  to  the  actuator  in  the  servo  must  be  constant. 
Tnio  means  that  tue  input  to  the, servo  amplifier  which  is  an  integrator  must  be  jjeroo 
In  order  for  this  to  be  true  must  equal  S  .which  in  turn  equals  <r  o.  In  other 
words  there  will  be  a  tracking  error  of  „  * 


Case  (2)  whs*1  the  missile  rotates  in  nitch  at  a  constant  rate  C  ijS  ),p  will 
rotatec  at  the  same  rate  but  in  the  opposite  direction  to. drive  S  (which  is  equal, 
to  @  +■  ft  )  to  zero.  In  .'the  absence  of  radomo  error  slope  this  can  be  accomplished 
resulting  in  no  tracking  error  being  required,  The  system  will  therefore  maintain 
both  £  *  and  S  at  zero  even  though  there  is  a  pitch  rotation. 


2  Guidance  Computer 


bo2ol  Guidance  Signal 


The  in  rnt  to  the  guidance  computer  should  be  .  It  has-been.  found 

in  the  previous  secti  ah  that  either  £  /  or  is  a  measure  of  in  the 

steady  state  conditions.  Qyna  icaliy,  however,  this  is  not  true..  The  transfer 
function  for  €l?  vrf.il  be  derived  first. 
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This  shows,  that  if'  £  were  used  fox*  guidance  cr  would  be  measured  ;&long 
with  cr  i  \V  and  ^  <*  There  will  also  be  a  time  lag* 

j  * 

The  expression  for  will  be  derived  next 
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This  again  contains  the  time  lag  and  a  pi ten  rate  coupling  term* 

If,  however,  the  combination  of  the  signals  g:  *HS>-  is  used  the  .following 
results: 

Differentiating  equation  l|o2~l 

*y  c jr(tr  ^l)  ( /*  ) *•  cr( r+0  ~  r  ft'  ( i*  ip*  )’*'*' ft 
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This  8fuia  result  was  obtained  in  section  2o)4*2  by  considering  only  the* 
definition  of  the  angles*  However  to  get  the  effect  of  reduced  £v  or  ,$•  the 
above  equations  must  be  used* 
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'  ♦  /  . 

The  ©-a  £  signal  contains  no  time  lags*  The  :  Heritable  radome  coupling 
tern  still  exists  however,. 

The  $'/*{■&  coupling  term  is  no  longer  present*  In  the  past  when  *  was 
used  for  guidance  the  y  Mk.  term  was  destabilizing#  l?or  the  DPN^k  system  K 
was  sss:  76  .  Analysis  showed  that  if  K  decreased  to  $0  there  will 

be  stability  problems,, 

•  1  i.  * 

With  the  &  signal  used  for  guidance  however  there  is  no  restriction 
on  the  value  of  K  dui*ing  guidanc  phases  of  the  flight*  The  maximum  value  for 
K  can  therefore  be  cbtermined  by  the  requirements  on  the  head  response  during 
the  nonguided  phased  of  flight*. 

p  /  *  / 

There  is  also  an  advantage  in  using  the  <$??£:  guidance  signal  from 
system  bias  consideration*  With  this  signal  only  the  bias  from  the  -rate  gyro 
appears  as  a  guidance  signal.  All  the  other  biases  such  as  df*  bias,  and  head 
servo  biases  result  in  a  constant  tracking  error  and  since  df  '  is',  differentiated 
the  constant  tracking  error  does  not,  appear  in  the  guidance  signal* 


lie 2* 2  Effect  of  'Radome  on  Guidance, 

•  /  «  j  ’ 

Using  O  as  the  g  'idance  signal  the  dosed  loop  transfer  function 
for  '°/dr  car.-  be  derived  assuming  a  constant  -/•?„ ,  (The  open  loop  transfer  function 
has'  already  Ween  derived  in  section  2*h*2), 
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The  steady  state  gain  of  the  system  has  been  modified  from 
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can  be  seen  from  the  equations  for  the  airframe  to  equal *  This 
means  that  the  larger  the  amount  of  o<  requiiad  to  produce  a  unit  y-  „ 
the  more  the  guidance  system  is  affected  by  raclome  error  slopes*  -or  a  given 
airframe  this  occurs  as  the  altitude  is  increased* 

The  actual  effect  of  f*  depends  on  its  sign*  If  it  is  negative  it 
might  cause  system  instability,  if  .it  is  positive,  the.  system  is  slowed  up  and 
also  causes  system  instability  if  the  autopilot  is  n^t  designed  properly* 

The  possible  means  of  coping  with  r  are  the  following: 

(1)  Discredit  large  r's  by  statistics* 

(2)  Artificially  bias  off  the  negative  r's  by  feeding,  into  the  autopilot 
a  term  porportional  to  f%ga  $  (explained  in  section  2Jj*2).  This 
results  in  a  situation  where  only  positive  t~  <s  need  to  be  dealt 
vdtho 

(3)  Design  a  compensating  network  to  stabilize  the  system  for  negative 
p's  *  A  lag-lead  network  has  been  found  to  be  feasible  when  the 
flight  condition  is  fixed* 

(U)  A  combination,  of  (2)  and  (3)* 

The  advantages  and  disauvantages  of  (2)  and  (3)  are: 

The  adv mintage  of  (2)  is  that  it  is  simple  and  the  resulting  system  can  be 
extremely  stable*  The  disadvantage  is  that  the  system  is  deliberately 
slowed  down* 

The  advantage  of  (3)  is  that  the  system  can  operate  with  both  positive  and 
negative  r's  and  has  r  «  0  as  the  nominal  condition*  The  disadvantage 
is  that  the  parameters  cf  the  lag-lead  netviork  will  probably  have  to  be 
varied  with  flight  conditions* 

The  Tartar  guidance  computer  was  designed  according  to  (2)*  The  autopilot 
was  designed  to  accomodate  a  (  r~  -h  hrM M  )  »  coB.  This  autopilot  along  with 

the  guidance  filter  of  _ J — TfiFrisrf)*'  »  howsvey  results  in  a  situs  c-Loa 

where  the  system  might  possibly  tolerate  negative  i~'s  ,  This  means  that 
2  T_f  *  Ta  <•  ?*  3  is  sufficiently  greater  than  ;VrT^  o  Or  in  the  Bode 
plot  of  Figure  9,  if  the  phase  is  changed  by  180°  the  system  is  still  stable  for 
certain  values  of  p  *  It  might  be  possible  therfore  to  use  something  less  than 
ft  feedback* 

^*2*3  Gui.dance  Filter 

The  guidance  filter  has  been  shown  to  be 
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In  order  to  shew  the  need  for  a  second'  order  filter,  the  total  transfer 
function  between  (f  noise  and  will  be  derived. 


therefore 
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In  order  for  f\>  to  remain  finite  when  &noise  is  white  noise, 

Gfr  (S)  must  be  at  least  secor.;  order*  The  ,?5  second  valueo  where  determined 
experimentaly*  These  values  can  be  subject  to  change  with  changes  in  the 
estimates  of  6.  noise  ex.ee  ct-?d» 
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IjoSoli  Effect  of  Head  Rate  Gyro-  Limitations, 

Thn  ^T0  nullvalue  appear?  directly  as  a  guidance  signal  bias. 

The  equivalent  Uc  due  to  this,  bias  can  be  computed  simply 

•-  nc  -  Sh 

The  g  sensitively  of  the  heac  rate  gyro  has  the  effect  of  introducing 
anotner  feedback  loop.  b 

The  guidance  signal  for  r  *  0  a  s  now 

4 Q~  b  €  b  q  n  ^  ^ 

•U& 

Khe»  K  is  the  g  sensitivity  o.f  the  head  rate  jy.  in  Vsec/g. 
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This  shows  t*,at  the  guidance  gain  and  all  the-  coefficients  of  the  guidance 

transfer  functions  are  modified  by  . . . L  -. _ _  b 

,  - 

AV,  ,  ,  ,svs' 

U  7Y7F  ^  1  aKd  K  “  °03  ^/sec/g  this  represents  12.% 
variation  on  the  gains  and  coefficients. 

iu2.5  Effect  of  Head  -Servo  Amplifier' Bias. 

a  bias  exists  at  the  input  to  the  servo  amplifier  there  is  no  affect 
on  th8  guidance.  However  during  the  head  positioning  phase  there  is  a  H  (£5) 
ieeaback  from  ^  to  servo  amplifier  in^at  and  the  effect  of  this  bias  is 

p  -  41  LA  -Hu>(si 

where  X\  is  the  unbalance  at  the  grid 

e~  1_ 

b  H  ( s) 
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■  This  shows  that  the  steady  state  error  due  to  J&L.  is  dependent  on 
the  zero  frequency  gain  of  H  (S)„ 


BLOCK  DIAGRAM  OF  TYPE  B  SYSTEM 
Fig,  2  (b) 
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APPENDIX  1* 

DERIVATION  OP  m®0J? 
FUHCTICS*  FOR  ELASTIC  BODY 


The  equation  for  an  elastic  Learn  is 

<n  i!  fXj< 

«fc'  <*x*  /£  ' ' 

Assume  the  solution  is  the  produet  of  &  funstios  of  ,y  and  *»v 

=;  <[&)/’(*) 

Substituting  into  equation  (1) 


F  ft)  d  "  £  w  t  (x)  xl  r 


dividing  through  by  f  {{)  /ly  fx  ^  ^ 

J...  JS-  £  iv)  i  a]  j/,.  <f  {%) 

'.'X 1  ~V'"" 


>1  (  A  ^  ^  (i) 

■  -  ciX 


JL~  JL~~  A  tf 
FkU  ci± 1 


Sloes  jp  f;o  does  net  vary  -sdth  -£  and  r-  (-:-)cio&8  not  vary  w».th  X  both 
aids.1  of  the  aquation  saust  be  equal  fce  a  constant 


**■  -1*  _  f~  (-'j 


-  T  W 


^  ,- 


-A.  JL  £  tr'  Sit)  ?!  i  ,1,1,  _  ,  »• 
MfiO  ✓x1 
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mltxplying  by  ■  and  jta&egtfat£ftg  fi'oa  0  to  L 


■  r'-  » 

y  ?*a)  ^  J*  -  j  %  I  /?£  »!«  «  ,'j)^ 


only  fcesa  tfiih  ^  xs.  &»  will  ysrayin 


/  /  /  *> 

/  $  £/>  -  /?.  /  j  <41 


/?•  -4  #  ^t***)^# 

i  '  at  ** 

■  fr;  »:*j  j  <y  D;)1  •//< 
«.  - 


The  nuiasrafcor  io  called  the  generalltied  and  ftho  tivcofelwita?-  fclio 

genstalissed  ss&ss 


liquation  {6}  hsc«sos 


)  /  -5/  /£)  d*  £(-5/  fa)  pv  t;!;- 

CTC  {  or  ^  '' 


?<■*  >  ^;/(lVl:/.i 

f  (.  *  »  ♦  '  4„-  > 


‘Hiere  aro  equations  oi*  the  fowi 


(8j  _ 4-,  ^  f.  ,  .J~  d.  lr,  u 

7n  i< )  p i  a)  <* «  *  ^ 4  f  {£')  2P- 


<•> .  £' 


>r  since  both  sides  of  tha  ;sque ttor  cof.  03  .•?«>•  'rcmal  to  a  constant.  ...  /.-y 


L  °  Xba  force 
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The  accelerometer  output  in  g*s  tharafora  bgcoaea 


hf/*)  -  5~;  M  j_(.^j 

(  S’ 3 ;-  j  ju<  s  32. 


The  rate  gyro  output  be coses 


f(dwy  -  sfe*)  L(^)  /w) 

:$*+?.  fu/ s  +  yn{  -«e* ) 


i-  for  tail  doflectior;  loading  only  bscoasa 


L  -  /v^/  Xs/i*Cf,ss 
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Subjests  Description  of  the  Coordinate  System  to  bo  Used  Haro  after  in.  the  Tartar 
Induced  Roll  Phenomena  Investigation 


IKTRODUCTIOH 

All  dynamics  study  of  the  Tartar  induced  roll  phenomena  sill  hereafter  us?  the 
coordinate  system  described  ia  this  memoo  This  system  differs  from  that  used 
heretofore  In  the  Tartar  studies,  free  that  used  by  APL  in  the  £>T£=>5  studies,  or 
that  usod  by  Conveir  personnel  in  the  STV=5  studies  o  It  is  balievcd,  ho  saver,  that 
this  is  the  system  that  has  bscn  agreed  to  be  everyone  concerned  as  the  one  to  usoo 

Only  those  features  of  the  system  that  are  pertinent  to  the  dynamics  studies 
t-o  be  performed  viU  be  presentedo  Those  studies  will  be  restricted  to  analyzing 
small  perturbations  in  yaw  and  roll  with  fixed  pitch  conditions o  The  equations  of 
motion  will  be  written  in  a  body  fixed  coordinate  system  The  orientation  of  this 
coordinate  system  in  the  body  will  depend  on  the  initial  roll  attitude  for  eny 
particular  phase  of  the  studya  The  aei’odynanic  forces  must  therefore  be  resolved 
to  that  body  fixed  coordinate  system  wlich  will  be  used  for  any  initial  roll  rtti~ 
tudeo  The  major  portion  of  this  memo  td32  be  involved  with  the  stability  derivatives 
required  for  these  aerodynamic  forceso  The  relationship  between  these  stability 
derivatives  and  those  for  the  forces  rosaLved  along  the  wind  tunnel  axis  is  shown*. 

The  equations  to  be  used  in  the  study  era  also  devolopedo 

DESCRIPTION  OP  COORDINATE  SSBTSi 

figure  1  shows  the  orientation  of  the  coordinate  system  in  the  missile  body<> 

The  origin  is  at  tho  misEUe  Co  go  The  sense  of  positive  rotations  are  also  shoisio 
These  conform  to  a  right  handad  systc^io  Also  shown  aro  the  orientation  of  tho  Z 
axis  with  respast  to  the  wings  for  the  cases  of  J$  «  0  and  0  -  •=h$°a 

figure  2  shows  a  positive  rotation  about  the  T  axis  from  XYZ  to  Tfl  €0  If  th© 
missile  velocity  vector  is  along  the  X  axis  aB  shown,  the  positive  rotation  produces 
a  positive  angle  of  attach  a*,  a  It  should  ba  noted  that  this  positive  <#?,  produces  a 
negative  force  along  the  Z  axis  (Fz)*  Also  with  this  sign  convention,  the  positive 
«*«  will  produce  a  negative  pitching  moment  (My)  for  a  staticly  stable  configuration™ 

Figure  3  chows  a  negative  rotation  about  the  Z  axis  fraa  XSZ  to  lC' Y^ZC  If  the 
missile  velocity  vector  is  along  tho  X  axic  us  shram,  tho  negative  I'irt&tion  produce® 
a  positive  side  slip  angle  (^?)«  This  positive  ^produces  a  negative  forct  along 
tho  Y  axis  (Fy)o  Positive also  produce  a  positive  yawing  moment  (Kg)  for  a 
staticly  stable  ccnfiguratioBo 

Tits  sign  convention  .for  yc  end used  ana  consistent  with  the  usual  dofinitiona- 
for  those  quantities 

(I)  sSja  o-x  -r~£ 


(;:)  sin  tjs 


Prepared  by 
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Figure  1}  shows  a  positive  rotation  (A  0)  about  the  X  axis  from  XI Z  to  X  ^ 
for  a  missile  pitched  at  au  initial  angle  of  attack  of  «:  0  and  zero  side  slip  angle* 

(This  would  correspond  to  the  conditions  under  which  the  wind  tunnel  tests  ero  run0) 
After  the  rotation^  <<  is  changed  and  thoro  is  also  a  resulting  ft  0 

The  magnitudes  of  <x  and  ft  cm  obtained  from  the  equations  (1)  and  (2)  defining 
and  ft  o  The  component  of  ?  along  the  t'  axis  (-*■/);  and  along  the  Z *  oxie  (gar) 
must  therefore  be  computed 0 

The  tranafozmation  matrix  for  a  rotation  about  the  X  axis  lag 

0  0  \ 

0  cos  0  sin  0 

0  ~sln  0  oos  0  y 

The  components  cf  V  before  the  rotation  were  as  follows 8 


U.  «  V  OOH  c<0 
0 

=  V  3ia*ry.0 

After  the  rotation  the  ccsnpoents  ara  as  follows  £ 
u!  CCS  <30  0 
'U*'  V  sin  ©o  o  sin  fi 
UrK  V  sin  «<  o  cfsa  0 

Using  tha  de.Cxttitican  for  szvi  (equations  3.  end  i) 
(3)  3lU  0<  i s  .'^£.  «  8*0  <a<(5)  cos  0 

(lj)  sin  «  sin  «3«r0  stai  0 

V 


In  tiie  dynanics  studies.?  t'w>  .  3  will  he  assaiaed  small,  and  0  wili  vary  a  small 
mount  &j#  fma  the  initial  orl  natation  of  the  Z  axis* 

The  equations  (3)  and  (U)  oeoca® 

(5)  <*<  a 

(6)  f  «  sia 
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It  should  be  noted  that  &C  and  fs  ara  in  the  planes  of  Z/6  V  and  Y^V* 
raspectivelyo 

figure  5  shows  the  sign  conventions  for  positive  deflections  of  tiio  tail 
surfaceso  In  the  0  »  0  roll  attitude  positive  &  produces  negative  Ffc  .and  negative 
My0  Positive  i?  produces  a  positive  Fy  and  a  negative  Kgo 

In  tbs  $  ®  °h$  roll  attitude  the  combination  of  <$f  and  ®l>  produces  a  positive 
Fy  and  a  nagativo  Hgo 

Figure  6  shows  the  differential  tail  deflections  ($' )  which  will  produce  a 
positive  rolling  moment  (Hjc)o  This  i3  defined  as  positive  $  0 

SlffiMAJCf  OF  COOIffllNATB  SYSTEM  DESCRIPTION 

The  information  presented  thuo  far  can  ba  susmarized  in  thD  following  table  0 
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cd  will  rosu&S'vfceewi® 

V 

mttmmwm  >■  i  . . 

„  ....  ,  .  -  f -  ..  _ 

<-*  & 

Produced  by  deflating  trailing  edge 

Produced  by  deflecting 

(deflections 

down 

trailing  edge  down  and  vi^x% 

dofinad 

Produces  -“Fg 

Produces  ^?2 

loolclng 

Produces 

Produces  ‘41y 

;  forward) 

Produces  m  Fy 

Produces  -Fy 

Produces  no 

1  Produces  <$$2 

I 

Produced  by  defi acting  trailing  edge 

Produced  by  deD.ectkng 

I  (deflections 

loft 

I  trailing  edge  down  and  left 

1  defined 

Produces  no  i:g 
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STABILITY  DERIVATIVES 

The  'tability  derivatives  reqci-ad  for  this  study  are  those  for  the  year  forcor 
yaw  dp-  ,  and  roll  momento  These  forces  and  moments  mist  be  referred  to  that 
bod'’  d  coordinate  system  which  is  used  for  any  particular  initial  roll  attitude 
(  Jon  aiders  do 

Figure  1  shews  the  body  fixed  coordinate  system  to  be  used  far  the  cas©  of 
fo  «  0  and  §b  =  *U5°o 

The  following  table  lists  the  stability  derivatives  and  tbs  signs  expectedo 
The  pitch  derivatives  are  included  though  they  are  not  needed  for  the  studies 
proposedo  The  signs  for  the  derivatives  can  be  deduced  from  the  informatics* 
presented  the  coordinate  system  descrlptiaa0  The  roll  yaw  coupling  derivatives 
with  respect  to  tail  surfaces  were  assumed  to  be  primarily  from  the  blanking  of  the 
upper  surfaces  when  the  missile  is  at  a  pnsitiva  «£  0» 


Table  of  Signs  Ejected  for  Stability  Derivatives 
(Assuming^  0  is  positive) 
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RATIONS  OP  HOTXOH 

The  required  equations  of  motion  for  this  study  ares 

o  . 

(7)  Fy  *  a  Or  <=*  pa  <•  ru) 

(8)  «  s*  Ig  o  (Iy  «  IX)  pq 

(?)  %»pix 

Making  the  following  substitutions 

p  =  0 

r  a  Y 

q  ®S<)  (constant  pitch  rats  associated  with  constant  angle  of  attack) 

o  * 

V  a  V  CO*  ft/& 

since  sin /& s  v/V  (equation  1) 
w  *  V  sin  <af  0  ( equation  2) 

u  *=  V  cos  of  0 
Equations  (7  9)  become 

(30)  Fy  «  a  (X  coa/j1^  ~  V  sia*r0  <>  cos  ©C0) 

(11)  %  «  Iz  #  ♦  (ly  -  IX)  /»  O0 

(12)  %  =  IXJ? 

>,  o  »  o 

whar®  Fy  is  in  lbcg  Hz  and  %{.  in  ftQ  Ibsg^  p  0e  y  „  ®0  are  in  rad/3*®* 
y>  ,  in  rad/ae®^ 

The  equations  fox’  tho  aerodyawnic  farces  and  acsacirts  aras 

(13)  Py  ^  3181 7\  S  I#  (Gy^4  ♦  Gay  iX  ♦  Cy^  $  ) 

Ul)  1%  «  3181  ?\  ed  K2  (Cj,/?  0  o  Cm  iy  *  Qn$  $  ) 

(15)  %  ~  1182. >.  od  (CX  $  <>  01(bp^  ^  ♦  Cx  .  ^  o  C^y  iy) 

6  ?V  2V  P 

The  iy  coefficients  are  defined  in  the  following  table o 


r  }  & » o 

4  «  45' 

i’X 

it 

if  -  i 

03U 

Gnf 

Cx£  -  Cyx. 

%X 

on,' 

cm;  “  °iu  • 

cli* 

°24?  ^  °li/ 

comxr. 

Porous 
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It  should  bo  noted  C(  for  the  0  =  case  equals  o70?  C(  uf  for  ths 
0  a  0  OBSOo 

Redefining  ths  sercdynsaic  coefficients  in  ths  following  manses’ 

A  *  3l|8l  A  S  H*  gel  Cy^s 

B  «  11*81  A  S  K2  £L&  Cm.’ 
nV 

C  a  11*81  A  8d  M2  ||a3  CH  ;9 
E  «  11*81  X  ed  H2  J&s3  Cgj,Y 

h 

F  «  11*81  A.  ad  M2  J&sl  *  C,. 

IX  5v  ^ 

'  27 

0«U*8XAedH2ikJ 

Ix  ft 

H«U*8lAsdi£  JKalCjL- 
IX  r 

X.  e  11*81  A  M  M2  Cli  Y 
«n*8iAs  a*  cY 

isV  0 

N  =  11*81 'AM  342  JLtJs 
*2 

Equations  (10  «  12)  bocoa« 

j  O  ■» 

(16)  ,S&I  Fy  «  ( A  o  «>  H  §  )  »•  co@^  sinenf  0  ^  Tcos  «<0 

BlV 

(17)  IZ&L&  0 

Iz  ^  *Z  $?«3 

(18)  2Z&22&  «  0  $  O  ffi  *  Rj$  ♦  kv  as  0* 

Ix  r 

For  these  equations  all  angles  ara  in  degress o 
The  control  equations  ero 

i-Yo  «*  %  (N  -  Hg)  «  Bg’JPo  fyy 
5©  55  -  (%  *  %  S  e  S2)  0 
where 


H  (gas)  ts 


V 


U.*fr%3  *H£) 


